arXiv: 1506.02613vl [astro-ph.SR] 8 Jun 2015 


Research in Astron. Astrophys. Vol.O (200x) No.O, 000-000 
http://www. raa-journal. org http://www. iop. org/journals/raa 


Research in 
Astronomy and 
Astrophysics 


Evolutionary Stages and Disk Properties of Young Stellar Objects 
in the Perseus Cloud 

Hong-Xin Zhang**^’^’^’^’^°, Yu Gao^’®, Min Fang®’^, Hai-Bo Yuan^, Yinghe Zhao®’®’^, 
Ruixiang Chang^, Xuejian Jiang®’®, Xiao-Wei Liu^, A-Li Luo^, Hongjun Ma®’®, Zhengyi 
Shao® and Xiaolong Wang®’® 

^ National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 
^ Department of Astronomy, Peking University, Beijing 100871, China; hongxin@pku.edu.cn 
^ Kevli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China 
Chinese Academy of Sciences South America Center for Astronomy, Camino El Observatorio 
#1515, Las Condes, Santiago, Chile 

® Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China; 
yiigao @ pmo. ac. cn 

® Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China 
^ Departamento de Lisica Teorica Universidad Autonoma de Madrid, 28049 Cantoblanco, Madrid, 
Spain; mfang.cn@gmail.com 

® Infrared Processing and Analysis Center, California Institute of Technology, MS 100-22, Pasadena, 
CA 91125, USA 

® Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, 
Chinese Academy of Sciences, Shanghai 200030, China 

Current address: Institute de Astroffsica, Lacultad de Lisica, Pontificia Universidad Catolica de 
Chile, Av. Vicuna Mackenna 4860, 7820436 Macul, Santiago, Chile 


Abstract We investigated the evolutionary stages and disk properties of 211 Young stel¬ 
lar objects (YSOs) across the Perseus cloud by modeling the broadband optical to mid- 
infrared (IR) spectral energy distribution (SED). Our optical gri photometry data were 
obtained from the recently finished Purple Mountain Observatory (PMO) Xuyi Schmidt 
Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC). About 81% of 
our sample fall into the Stage II phase which is characterized by having optically thick 
disks, and 14% into the Stage I phase characterized by having significant infalling en¬ 
velopes, and the remaining 5% into the Stage III phase characterized by having optically 
thin disks. The median stellar age and mass of the Perseus YSOs are 3.1 Myr and 0.3 Mq 
respectively. By exploring the relationships among the turnoff wave bands Atumoff (long- 
ward of which significant IR excesses above the stellar photosphere are observed), the 
excess spectral index ctexcess as determined for A > Atumoff, and the disk inner radius Rin 
(determined from SED modeling) for YSOs of different evolutionary stages, we found 
that the median and standard deviation of aexcess of the YSOs with optically thick disks 
tend to increase with Atumoff, especially at Atumoff > 5.8/im, whereas the median frac¬ 
tional dust luminosities Ldust/L* tend to decrease with increasing Atumoff - This points to 
an inside-out disk clearing of small dust grains. Moreover, a positive correlation between 
aexcess and i?in was found at aexcess ^ 0 and Rin ^ 10 x the dust sublimation radius 
Rsub, irrespective of Atumoff, Ldust/L* and disk flaring. This suggests that the outer disk 
flaring either does not evolve synchronously with the inside-out disk clearing of small dust 
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grains or has little appreciable influence on the spectral slopes at A ^ 24/im. About 23% 
of our YSO disks are classified as transitional disks, which have Atumoff > 5.8pm and 
Ldust /L* > 10“^. The transitional disks and full disks occupy distinctly different regions 
on the Ldust/L* vs. aexcess diagram. Taking Ldust / L* as an approximate discriminator of 
disks with (>0.1) and without (<0.1) considerable accretion activity, we found that 65% 
and 35% of the transitional disks may be consistent with being dominantly cleared by 
photoevaporation and dynamical interaction with giant planets respectively. None of our 
transitional disks have aexcess (<0.0) and Ldust/L,t (>0.1) values that would otherwise 
be suggestive of disk clearing dominantly by grain growth. 

Key words: stars; formation - stars; low-mass - stars; pre-main sequence - individual; 
Perseus Cloud - circumstellar matter - protoplanetary 


1 INTRODUCTION 

The formation and early evolution of stars are among the central problems in Astrophysics. Young stellar 
objects (YSOs), which are primarily identified as pre-main-sequence (PMS) stars by the presence of 
infrared (IR) excess arising from circumstellar disks or surrounding envelopes (e.g. Allen et al. 2004; 
Greene et al. 1994; Lada 1987), have been extensively studied in nearby star-forming regions (e.g. 
Taurus; Luhman et al. 2010; NGC 1333; Winston et al. 2010; IC 348; Muench et al. 2007; a Ori; 
Hernandez et al. 2007; Tr 37; Sicilia-Aguilar et al. 2006; NGC2362; Currie et al. 2009; Lynds 1630N, 
1641; Fang et al. 2009,2013). Studying the circumstellar environment, either disks or envelopes, around 
YSOs of different mass is essential to understanding the formation of both stars and their planetary 
systems. 

YSOs are traditionally categorized into four classes or evolutionary stages based on the spectral 
index a (dlog(AF(A))/d log (A)) of their near- to mid-IR spectral energy distributions (SEDs; e.g. Andre 
et al. 1993; Greene et al. 1994; Lada 1987). The youngest Class 0 objects are only visible in far-IR to 
submm wavelengths, and they are thought to have envelope mass that exceeds the central stellar mass; 
Class I YSOs (a > 0.3) are characterized by rising mid-IR SEDs, and may be still in an envelope 
collapse stage but have the central stellar mass exceeding the envelope mass; Class II YSOs (—1.6 < 
a < —0.3) have SEDs that peak at near-IR wavelengths, decrease at longer wavelengths that is much 
more gradual than what expected for a stellar photosphere, and they agree well with PMS stars with 
circumstellar accretion disks; Class III YSOs (a < —1.6) have little or no IR excess, and are thought 
to be in the disk dissipation stage with very little or no circumstellar material. In addition, Greene et al. 
(1994) introduced an additional pLAT-spectrum class (—0.3 < a < 0.3) which has spectral indices in 
between Classes I and IT 

The star formation process is generally accompanied by the formation, evolution and dispersal of 
circumstellar protoplanetary disks which are believed to be the sites of planet formation. In particular, 
optically thick full disks are usually found in the Class II YSOs, whereas the evolved or anemic optically 
thin disks are usually identified with the Class III YSOs. A lot of important information about the evo¬ 
lutionary stages and disk properties of YSOs is encoded in the multi-wavelength SEDs (e.g. Robitaille 
et al. 2006, hereafter R06). Eor instance, the optical to near-IR bands offer important constraints on the 
properties of the central source, such as the temperature and bolometric luminosity; the near- to mid-IR 
bands provide a crucial constraint on the inner (from a few AU to tens of AU) disk properties; the far-IR 
to submm bands give strong constraints on the mass of disks and envelopes (e.g. Andrews & Williams 
2005). 

As currently the most active site of low- to intermediate-mass star formation within ^300 pc of 
the Sun, the Perseus cloud (M ~ 4.8x10^ Mq; Evans et al. 2009) region is an ideal laboratory for 
studying the formation and early evolution of low- to intermediate-mass stars (e.g. Bally et al. 2008) and 
the circumstellar disks. Recently, the Spitzer telescope observations, especially through the “Cores to 
Disks” legacy project (c2d; Evans et al. 2003), have led to the identification of over 400 YSOs (mostly 
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Classes I and II) toward the Perseus cloud. In addition, systematic submm continuum surveys of the 
Perseus region with SCUBA (Hatchell et al. 2005) and Bolocam/CSO (Enoch et al. 2006) have led 
to the confirmation of over 100 protostellar or starless submm cores, and about one-third (one-fifth) 
of these cores were classified as Class 0 (Class I) YSOs. About two-thirds of the Perseus YSOs are 
associated with the two major young clusters NGC 1333 and IC 348, and the remaining YSOs are either 
associated with other much smaller clouds, such as Barnard 5, Barnard 1, L1455 and L1448, or sparsely 
distributed across the whole Perseus cloud region (e.g. Evans et al. 2009; J(/)rgensen et al. 2007). 

A systematic investigation of the evolutionary stages and disk properties of the Perseus YSOs 
with the optical-to-IR SEDs is still lacking. Moderately deep broadband gri imaging data were re¬ 
cently obtained through PMO’s Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center 
(XSTPS-GAC; Liu et al. 2014; Zhang et al. 2013, 2014; Yuan et al. 2015, in preparation). In this paper, 
we combined the gri data with the IR data from 2MASS, Spitzer, and WISE in order to study the 
physical properties of the central stellar sources, the evolutionary stages and inner disk properties of the 
Perseus YSOs. Euture spectroscopic data from LAMOST on most YSOs will further detail and enhance 
the broad-band characterization offered in this paper. Section |2] introduces the data and YSO sample 
analyzed in this work. The color-magnitude diagrams are presented in Section[3 Section |4]presents the 
results from SED modeling, such as the central stellar masses, ages, and the evolutionary stages of the 
YSOs. An investigation of the excess dust emission and disk geometry parameters, such as the disk 
inner radii and outer disk flaring, and implications on the dominant disk clearing processes are given in 
Section|5] A brief summary of the main results in this work is given in Section|6] 

2 SAMPLE AND DATA 

2.1 Parent Sample of Perseus YSOs 

The most recent census of Perseus YSOs was done by Hsieh & Lai (2013, HL13) using the photometric 
data from Spitzer c2d legacy project (Evans et al. 2009), which carried out a wide-field imaging survey 
of five nearby low-mass star-forming clouds (Serpens, Persues, Ophiuchus, Lupus, and Chamaeleon) 
with both IRAC and MIPS instruments onboard Spitzer. Instead of simply relying on a cut on one or 
two color-color and color-magnitude diagrams to separate YSOs from main-sequence stars and back¬ 
ground galaxies, HL13 identified YSOs in a multi-dimensional magnitude space. In particular, HL13 
used data from the Spitzer SWIRE survey of the ELAIS N1 extragalactic held (Surace et al. 2004) to 
acquire a control sample for background galaxies, and this control sample was used to define the regions 
occupied by galaxies in the multi-dimensional magnitude space. The readers are referred to HE 13 for 
more details about the YSO identification procedure. 

In total, HL13 identified 469 Perseus YSOs over 3.86 deg^ covered by the c2d survey. Adopting a 
distance of 250 pc for the Perseus cloud, 3.86 deg^ corresponds to about 73.6 pc^ (Evans et al. 2009). 
Among the 469 YSOs, 21% was classified as Class O/I sources, 10% as Class Elat sources, 58% as 
Class II sources, and 10% as Class III sources based on the 2MASS Kg to MIPS 24/rm spectral indices 
a. We note that 429 of the 469 YSOs have detections in at least 3 IR bands, and thus the identification of 
these 429 YSOs in the multi-magnitude space should be more reliable than the rest 40. In the following, 
the 429 YSOs will be regarded as the parent sample, and our subsample selection and analysis will be 
based on these 429 YSOs. 

2.2 Data 

2.2.7 Broadband gri Photometry from XSTPS-GAC 

Prom the fall of 2009 to the spring of 2011, the XSTPS-GAC observing project carried out an imaging 
survey toward the Galactic anti-center in SDSS gri bands with the PMO’s Xuyi 1.04/1.20m Schmidt 
Telescope. The survey covers the sky area from RA ^ 45° to 135° and DEC ^ —10° to 60°, plus 
an extension of ^ 900 deg^ toward the M31/M33 direction. With an exposure time of 90 seconds, 
the survey reaches rum ^ 19 in the r band at lOcr for point sources. The astrometry (accurate to ^ 


4 


H.-X. Zhang, Y. Gao, & M. Fang et al. 


0.1") was calibrated against the PPMXL catalog (Roeser et al. 2010), and the PSF-fitting photometry 
was calibrated against the SDSS DR8 catalog using the overlapping sky area with an accuracy of 2%. 
Given the importance of optical bands in constraining the properties of central stellar sources of YSOs, 
XSTPS-GAC point sources with S/N > 2 (riim ~ 21 mag) will be used in this work. 

2.2.2 Spitzer Data from the c2d Project 

As mentioned above, the Perseus cloud has been observed by the c2d project in the Spitzer IRAC 3.6 
(IRl), 4.5 (1R2), 5.8 (1R3), and 7.9 (1R4) pm and MIPS 24 (Ml), 70 (M2), and 160 (M3) pm bands. 
All data, including imagery and point-source photometry (through PSF fitting) catalogs for IRAC, Ml 
and M2 were processed and released by the c2d team. In this work, we used the HREL (high reliability) 
source catalog provided by the c2d projecQ. 

2.2.3 2MASS and WISE Data 

The JHKs photometry was taken from the 2MASS Point Source Catalog which reaches a Kg-hwd 
limiting magnitude of 14.3 mag at lOcr. The Wide-held Survey Explorer (WISE) survey mapped the 
whole sky in four IR broad bands, i.e. 3.4 (Wl), 4.6(W2), 12(W3) and 22(W4) pm, with a Sa limit¬ 
ing magnitude of 16.6, 15.6, 11.3 and 8.0 mag respectively for the four bands. In this work, we used 
the ALLWISE Source Catalog which includes enhanced photometric sensitivity and accuracy, and im¬ 
proved astrometric precision compared to the 2012 WISE All-Sky Data Release. 

2.3 Our Working Sample 

In this work, we selected a subsample of 211 Perseus YSOs from the HL13 parent sample. The 211 
YSOs were selected by cross-matching the HL13 catalog with all the above data sets, with a requirement 
that each source should have JHKs, IRAC or WISE, Ml or W4, and at least one optical band available. 
Among the 211 YSOs, 102 have p-band detections with S/N >2, 151 have 7’-band detections, and 198 
have i-band detections. We point out that 78% (99%) of the p-band detections have S/N > 10 (5), 85% 
(99%) of the r-band detections have S/N >10 (5), and 94% (100%) of the /-band detections have S/N 
> 10 (5). In addition, 27 of our sample YSOs have M2 detections. Optical photometry of the 211 YSOs 
is given in Table [1] 

Spatial distribution of the Perseus YSOs is shown in Figure [T] the 110 GHz ^^CO integrated in¬ 
tensity map from the The Coordinated Molecular Probe Line Extinction Thermal Emission Survey of 
Star Forming Regions (COMPLETE, Goodman et al. 2005; Ridge et al. 2006) project is overlaid for 
comparison. As already known, most the Perseus YSOs are associated with the two major clusters IC 
348 and NGC 1333. In particular, 83 of the 211 YSOs are within 15' (^ 1.7 pc at a distance of 320 pc, 
e.g. Belikov et al. 2002; Evans et al. 2009; Strom et al. 1974; de Zeeuw et al. 1999) of IC 348, and 43 
are within 15' (^ 1.3 pc at a distance of 250 pc, e.g. Evans et al. 2009) of NGC 1333. 

Eigure|2]shows histograms of the IRl mag and the spectral indices a(iTs—Ml) for the parent sample 
and our working subsample. As can be seen, our 211 YSOs are expected to be statistically unbiased at 
least at IRl ^ 10 mag, which would correspond to a stellar mass of ~ 0.9 Mq at an age of ~ 3 Myr 
for a distance modulus (to — M)q = 7.5 (corresponds to a distance of 320 pc for IC 348), according to 
the PMS evolutionary tracks of Baraffe et al. (1998). Moreover, since hot dust of the circumstellar disks 
may contribute significantly to the IRl emission, our subsample of YSOs may be statistically unbiased 
down to M* of slightly below 0.9 Mq. In addition, most of our YSOs have a ^ 0.0, implying that our 
subsample is dominated by Classes Flat, II and III YSOs. The spectral index a(iTs— Ml) (Evans et al. 
2009), which quantifies the spectral slope from Ks to Spitzer 24/rm, was obtained from a linear fit to 
logarithms of all available photometry between Kg and Ml. Note that for sources without Ml data we 
used the W4 data instead for determining the spectral indices. 

* http://irsa.ipac.caltech.edu/data/SPITZER/C2D/ 

^ http://wise2.ipac.caltech.edu/docs/release/allwise/ 
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Fig. 1 Spatial distribution of the Perseus YSOs is over plotted on the FCRAO 110 GHz 
integrated intensity map (greyscale, FWHM ~ 46”) from the COMPLETE project. The small 
red circles mark the 211 YSOs studied in this work, and the blue crosses mark the parent 
sample of 429 YSOs which have at least 3 band detections in the IR wavelength range from 
2MASS J to MIPS 24pm. Several well-studied clusters or cores are also annotated in the 
figure. 



Fig. 2 Histograms of the IRAC 3.6pm (IRl) magnitude (left) and spectral index a (ATg—Ml) 
(right) of wavelength ranges from Kg to MIPS 24pm for the parent sample (open) and our 
working subsample (hatched) of YSOs. 
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Fig. 3 J — Ks vs. Kg (left) and / — J vs. J {right) color-magnitude diagrams. The left 
panel shows the distribution for all of the 211 YSOs (filled circles) studied in this work, and 
the right panel shows the distribution for 198 YSOs with i-band detection. Overplotted are 
stellar evolutionary tracks of Baraffe et al. (1998) for a stellar mass range of 0.02 — 1.4 Mq 
at three different ages (1,3, and 30 Myr). The black arrow in each panel marks the 5-mag 
visual extinction vector, assuming the Fitzpatrick (1999) extinction law with Ry = 3.1. The 
evolutionary tracks and extinction vectors shown in the right panel are the same as that in the 
left panel. A distance modulus (m — M)o = 7.0 for the Perseus YSOs is adopted. 


3 COLOR-MAGNITUDE DIAGRAMS 

The color-magnitude diagrams for our sample are shown in Figure |3] The evolutionary models for low- 
mass stars and brown dwarfs from Baraffe, Chabrier, & Allard (1998) are also plotted in Figure|3]to be 
compared with our data. When plotting the evolutional models in Figure |3 the JHK photometry on 
CIT system as provided by Baraffe et al. (1998) was transformed to the 2MASS photometric system, and 
our SDSS i magnitude was transformed to the Cousins I magnitude using the transformation equation 
determined by Lupton (2005). The transformation equation of Lupton (2005) involves Cousins I, SDSS 
r, and SDSS i. Among the 198 YSOs that have i-band detection, 60 do not have ?’-band detection. To 
put these 60 YSOs on the I — J vs. J diagram (right panel of Figure|3ll, we adopted a median I—i = 
—0.82, as determined from the YSOs with both r and i detections, to transform SDSS i to Cousins I. 

In Figure [3 most of our YSOs are redder than the pure stellar photosphere emission. A recent 
study by Chen et al. (2015) found a mean visual interstellar extinction of 1 mag toward the Perseus 
region, which is insufficient to explain the red colors of most YSOs, especially for their distribution 
on the J — Kg vs. Kg diagram. Therefore, as expected, hot dust emission from the inner circumstellar 
disks of YSOs contributes significantly to the Kg band. The comparison with theoretical evolutionary 
tracks implies that the masses of our YSOs are mostly above the substellar limit (^ 0.08 Mq). The fact 
that evolutionary tracks at different masses and ages are well separated on the color-magnitude diagram 
involving / band data suggests the importance of optical bands in constraining the properties of the 
central stellar sources of YSOs. 
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4 SED MODELING 

4.1 The Method 

With the broadband SEDs in hand, we used the online SED fitting tool developed by R06 and Robitaille 
et al. (2007) to extract the relevant physical properties of YSOs and their circumstellar disks. This online 
fitting tool offers the possibility of fitting YSO SEDs with a precomputed grid of 200,000 synthetic 
SEDs computed at 10 viewing angles. The model SEDs account for the contribution from central stellar 
photosphere emission, circumstellar disks, and infalling envelopes. In particular, the stellar photosphere 
emission is modeled with two parameters, i.e. stellar luminosity and temperature; The disk is treated 
as a standard flared accretion disk and the resultant emission is modeled with six parameters, i.e. the 
disk mass 0.001 - 0.1 Mq), inner radius, outer radius (1 - 10,000 AU), accretion rate, scale height 
factor and flaring angle; The envelope emission is modeled with four parameters, i.e. envelope accretion 
rate, outer radius, cavity density (10“^^ - 8 x 10“^° g cm“^) and cavity opening angle. In addition, the 
central stellar masses (0.1 - 50 Mq) and ages (0.001 - 10 Myr) were constrained by comparing the 
stellar luminosity and temperature with the PMS evolutionary tracks of Bernasconi & Maeder (1996) 
and Siess et al. (2000). 

Before proceeding to the SED modeling for our data, we point out some limitations of the R06 SED 
models (Robitaille 2008) that may be relevant to our current work. Eirstly, the models do not include 
the case for multiple central stellar sources, which can affect the size of the disk/envelope inner holes 
and thus influence the near- to mid-IR emission. Secondly, there exist several sets of different PMS 
evolutionary tracks in the literature, besides the Siess et al. tracks as adopted by the R06 SED models, 
other popular PMS tracks include Swenson et al. (1994), D’Antona & Mazzitelli (1997), Baraffe et al. 
(1998), Palla & Stabler (1999),Yi et al. (2003), and Dotter et al. (2008). Adopting different tracks can 
lead to systematic differences on the fitted stellar parameters (e.g. Pang et al. 2013; Hillenbrand et al. 
2008), and the systematic effects are especially significant for sub-solar mass stars at young ages. In 
particular, uncertainties of age estimation from different tracks for sub-solar mass stars can be up to 
0.75 dex at young ages (< 10 Myr, Hillenbrand et al. 2008). Thirdly, the dust opacity law assumed in 
the models may not be accurate, which would affect the determination of disk/envelope accretion and 
mass. 

When fitting the SEDs, a 5% absolute flux calibration uncertainty was added in quadrature to the 
gri uncertainties, a 10% uncertainty was added to JHKg and IRAC data uncertainties, and a 20% 
uncertainty was added to the M2 data uncertainties (Evans et al. 2009). In addition, when both IRl, IR2 
and Wl, W2 data are available, IRl and IR2 were used in the fitting due to the higher resolution of 
IRAC data. An aperture of 10" was used in the fitting. In addition, the distance to YSOs was allowed to 
vary from 0.2 to 0.35 kpc, and the foreground interstellar extinction Ay was allowed to vary from 0.3 
to 30 mag, with the lower limit of Ay being chosen based on the Perseus extinction map as determined 
by Chen et al. (2014). Besides the best-fitting model parameters, all the subsequent well-fit models 
with reduced x? — Xr best < 2 were used to define the minimum and maximum acceptable physical 
parameters. 

4.2 The Results 

The range of wavelength coverage determines what physical parameters can be constrained from SED 
modeling. A thorough investigation about how the wavelength range of data affects the determination of 
different physical properties of YSOs was given by R06. Given our wavelength coverage from optical 
to MIPS 24/im (or WISE 22/im), we expect to roughly constrain the central stellar source luminosity, 
extinction, and the circumstellar disk luminosity. Although subject to much larger uncertainties than 
constraints from spectroscopic data, the central stellar masses and ages can still be roughly constrained 
from broadband SED modeling to statistically investigate a large sample, like the one presented in this 
work. Moreover, while the masses of the circumstellar disks and envelopes cannot be reliably con¬ 
strained unless one have far-IR to submm data, SED modeling for wavelength ranges short of far-IR 
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can still be used to statistically constrain the evolutionary stages of YSOs. R06 found that at least 3 
different evolutionary stages of YSOs can be statistically distinguished based on the fitted stellar masses 
normalized envelope accretion rates (Menv/-^*) and disk masses (Mdisk/-^*)- In particular, the Stage 

I YSOs have significant infalling envelopes and are defined by having Menv/-^* > 10“® yr“^; Stage 

II YSOs have optically thick disks and are defined by having M^nv/M^ < 10“® yr“^ and Mdisk/-^* > 

10“®; Stage III YSOs have optically thin disks and are defined by having both < 10“® yr“^ 

and Mdisk/-M* < 10“®. Lastly, the near- to mid-lR SEDs are also sensitive to the disk properties, such 
as the disk inner radius and disk flaring. 

Figure 0] shows the SEDs of the 27 YSOs which have M2-band detections and at the same time at 
least one optical band available. The black solid curve in each panel of Figure |4] is the best-fit model 
SED, and the grey solid curves represent all subsequent well-fit models with reduced X?~X?best 
In addition, SEDs of the best-fit stellar photosphere emission (corrected for both the interstellar and 
circumstellar extinction) are overplotted as dashed curves. By calculating likelihood estimator e“^r/^ 
for each well-fit model with Xr ~ X? best < 2 for a given YSO, we construct the probability density 
function (PDF) and the correspondent cumulative distribution function (CDF) for parameters such as 
stellar masses M*, ages and disk inner radius i?in. The most probable value for each parameter refers to 
the median of the corresponding PDF, and the confidence interval is defined as covering the central 95% 
of the CDF. In what follows in this section, we will present the results for M*, ages, and the evolutionary 
stages as identified based on the stellar masses normalized disk masses and envelope accretion rates. 
Discussion about the disk geometry parameters from SED modeling and fractional dust luminosity 
Fdust/F*, where Ldust (in units of Lq) is equal to integral of the best-fit stellar photosphere subtracted 
SEDs, will be presented in the next section. SED modeling results for some relevant parameters, such 
as M*, ages, and i?in, are listed in Table|2 

4.2.1 Stellar Mass Distribution of the Central Stellar Sources 

Stellar mass histogram of our sample is shown in the left panel of Figure |5] As pointed out above, our 
sample is expected to be statistically unbiased at M* ^ 0.9 Mq. We overplot the Salpeter stellar initial 
mass function (IMF; Salpeter 1955) which was scaled to have the same number of stars at M* > 0.9 Mq 
with our YSO sample. The error bars in the histogram represent the Poisson noise from number counts. 
It can be seen that the mass distribution of our YSOs at M* ^ 1 Mq is consistent with the Salpeter IMF 
within the uncertainties. Note that an extended star formation history for the Perseus region might make 
it not straightforward to compare the accumulated present-day mass function with the simple Salpeter 
IMF. Although our sample may be subjected to significant incompleteness bias below 1 Mq, we note 
that a flat and broad mass distribution from sub-solar to sub-stellar mass limit, as found in our sample, 
is in general agreement with previous studies of low-mass clusters such as IC 348 (e.g. Luhman et al. 
2003a; Muench et al. 2003), NGC 1333 (e.g. Wilking et al. 2004; Greissl et al. 2007),Trapezium (e.g. 
Muench et al. 2002) and other nearby clusters (e.g. Andersen et al. 2008; Hillenbrand & Carpenter 2000; 
Luhman et al. 2000; Lucas et al. 2005; Luhman 2007; Levine et al. 2006; Moraux et al. 2003; Slesnick 
et al. 2004; Scholz et al. 2009; Weights et al. 2009). The median M* of our YSOs is ~ 0.3 Mq. 

Stellar mass distributions of YSOs within 15' radius of each of the two major clusters IC 348 and 
NGC 1333 are shown in the left panel of Figure |6] The median stellar masses of YSOs in IC 348 and 
NGC 1333 are-0.3 Mq. 

4.2.2 Age Distribution of the Central Stellar Sources 

The age histogram for the whole sample is shown in the right panel of Figure 0 and the age histograms 
for each of the two major clusters (again defined with a 15' radius) are shown in the right panel of Figure 
|6] The median stellar age of the whole sample is — 3.1 Myr, and the median age for YSOs in IC 348 
and NGC 1333 is — 2.8 and 2.5 Myr respectively. A relatively younger age of NGC 1333 than IC 348 is 


(ergs/cm^/s) XF^ (ergs/cm^s) (erga/cm^/s) XF^ (erga/cm^/s) XF^ (ergs/cm^/a) XF^ (ergs/cm^/a) XF^ (erga/cm^/s) 


YSOs in the Perseus Cloud 


9 


■ c2dJ032741.5+302017 

~E 

pio"” 

c2dJ032800.1+300347 

~E 10-’“ 

lio-’’ 

c2dJ032834.5+310051 

» i.m)' 


c2dJ032B50.6+304245 

X (/^m) 


X Oim) 




X (Mm) 

C2gj032929.3+311835 


c2dJ032954.0+312053 


f' c2dJ033001.9+303529 


c2dJ033035.9+303024 

> (»m) 


X (Mm) 




X (Mm) 

- C2gj033037.0+30312S 

/ m,, . 


c2dJ033044.0+303247 


c2dJ033052.5+305418 


C20J033110.7+304941 

X («m) 


X (Mm) 


X (Mm) 


X (Mm) 

c2dJ033114.7+304955 

\ 

~E 

c2dJ03311B.3+304940 

~E ,p-,o 

c2dJ033234.0+310056 

~E ,(,-10 

c2dJ033312.8+312124 

X (jim) 


X (Mm) 


X (Mm) 


X (Mm) 

c2dJ033330.4+311051 

~E 

10"" 

c2dJ033341.3+311341 

i 

10-'“ 

c2dJ033401.7+311440 

, 0-10 

c2dJ033449.8+31 1550 

X (jim) 


X (Mm) 


X (Mm) 


X (Mm) 

c2dJ034109.1+31443S 

"^E ,g-,, 

C20J034155.7+314B1 


c2dJ034157.4+314B37 

~£ 

c2dJ034219.3+314327 

X (urn) 


X (Mm) 


X (Mm) 



-• - c2dJ034344.5+314309 

■- ' / ''' \ 

~E 10"'“ 

i' 

c2dJ034516.3+320620 

~E 

P 

c2dJ034520.5+320634 



. . 1. 1 



lU 





X X Oim) X Oim) 


Fig. 4 SEDs of 27 Perseus YSOs. Among our whole sample, these YSOs have at least one 
optical band, JHKg, IRAC or WISE, MIPS 24/rm or WISE 22/j,m, and MIPS 70/j,m avail¬ 
able (black points). The black solid curve in each panel represents the best-fit model SED 
of Robitaille et al. (2007), and the grey curves represent all subsequent well-fit models with 
Xr “ Xt best 2. The dashed lines illustrate the SEDs of stellar photosphere in the best-fit 
model, as it would appear to be in the absence of circumstellar dust. 
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Fig. 5 Histograms of the stellar masses (left) and ages (right) for the whole sample of YSOs. 
The filled red circles in the left panel represents the Salpeter IMF that is scaled to have the 
same number of observed stars more massive than 0.9 Mq. 



M^(Mq,) Age(Myr) 

Fig. 6 Histograms of the stellar masses (left) and ages (right) for the two major clusters IC 
348 (thick red) and NGC 1333 (thin blue). YSOs within 15' radius of each of the two clusters 
are regarded to be associated with the cluster. The filled red (blue) circles in the left panel 
represent the Salpeter IMF for IC 348 (NGC 1333) that is scaled to have the same number of 
observed stars more massive than 0.9 Mq. 
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Fig. 7 Classification of evolutionary stages of our sample based on the fractional envelope 
accretion rates and disk masses. 


in line with previous studies, and our age estimate is consistent with previous studies of YSOs in these 
two clusters (e.g. Herbig 1998; Luhman et al. 2003b; Lada et al. 2006; Winston et al. 2009). 

4.2.3 Uncertainties on Stellar Parameters from SED Modeling 

Determination of the masses and ages of central stellar sources relies on a reasonably accurate con¬ 
straint on the effective temperature Te®. While it is reasonable to statistically explore the distribution 
of masses and ages determined from broadband SED modeling for a large sample, results for individ¬ 
ual sources may be subject to large uncertainties. In principle, T^s can be accurately constrained by 
photospheric absorption lines from optical or near-IR spectroscopy. By comparing our SED-based and 
the spectroscopy-based TeS for 75 IC 348 YSOs that have spectroscopic observations in the literature 
(e.g. Luhman et al. 2003b; Lada et al. 2006; Muzerolle et al. 2006; Muench et al. 2007), we found 
a median and standard deviation of Tefi^ED—Tea,spec of 71 (~ 2%) and 257 K 7%) respectively 
for the 35 objects with Ay < 4 mag and Tea spec < 5000 K, and a median and standard deviation of 
Tea, SED— Tea, Spec of 68 2%) and 434 K 12%) respectively for the 31 objects with Ay > 4 mag 

and Teff.Spec < 5000 K. In addition, the remaining 9 objects with Tea,Spec > 5000 K have a median 
and standard deviation of Tea, sed—T ea,Spec of —1333 (~ 24%) and 1229 K (~ 16%) respectively. 
According to the theoretical evolutionary tracks of Baraffe et al. (1998), for a PMS star with Tea of 
3336 K and an age of 3 Myr, which corresponds to a stellar mass of 0.3 Mq, an overestimation of Tea 
by ~ 350 K (~ 2-3 subclasses in spectral type) at a given luminosity can result in an overestimation of 
age and mass by factors of 3 and 2 respectively. 

4.2.4 Evolutionary Stages of YSOs 

Similar to Povich et al. (2013), for every YSO, we calculated the accumulated probability 
(TstageOcX^modei i 6“^*''/^) of it being in each of the three Stages (i.e. Tstagei, Tstageii, Tstageiii) based 
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Fig. 8 Breakdown of different evolutionary stages into different q;(Xs—M l)-based Classes. 


on all the well-fit models with Xr ~ X?best < 2. A YSO is uniquely classified as a Stage I, II, or III 
object if the normalized Pstage > 0.67. The result of our classification is presented in Figure]?] There are 
5% YSOs that can not be classified as either Stage I, II, or III if adopting the 0.67 probability threshold. 
These 5% objects were classified as evolutionary stages that have the highest accumulated probability. 
The classification of some YSOs into Stage I phase may be subject to relatively large uncertainties. 
This is because the wavelength coverage of our SEDs is mostly limited to ^ 24/im, shortward of which 
the contribution of excess emission from disks dominates over that from the cool infalling envelopes. 
Moreover, we note that some of our Stage I YSOs with low IR excess luminosities may be genuine 
Stage II YSOs with edge-on optically thick disks. 

As can be seen from Figure|7| our sample is dominated by Stage II YSOs. Moreover, the fractions 
of YSOs in different stages are similar for IC 348, NGC 1333, and the other regions. As mentioned in 
the Introduction section, YSOs have been historically grouped into three or four classes based on the 
spectral index a determined over the wavelength range from ^ 2 to 20 pm. YSOs of different classes 
are thought to be in different evolutionary stages (see above for references). R06 showed that there is a 
general correspondence between the modeling-based “Stages” and a-based “Classes”, in the sense that 
Stage I is expected to include the Class O/I, Stage II is analogous to Class II, and Stage III to Class III. 
However, as a set of purely empirical criteria, the Class scheme can be sometimes misleading. 

Figure [Sjpresents the breakdown of each Stage into different Classes. As is shown, a vast majority 
(94%) of Class II objects are grouped into the Stage II, and the majority (85%) of Class I objects are 
grouped into the Stage I. It is noteworthy that the dominant physical Stages for Class Flat YSOs are 
uncertain, with about 43% being in Stage I and the remaining 57% in Stage II. Fikewise, the dominant 
physical Stages for Class III YSOs are also uncertain, with about 46% of them being in Stage II and the 
remaining 54% in Stage III. 
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Fig. 9 Distribution of ctexcess vs. the wavelength bands Atumoff longward of which IR ex¬ 
cesses are observed. Stages I, II and III YSOs are shown separately in the left, middle and 
right panels. Median of aexcess at each individual Atumoff of different Stages was shown as 
red triangles. Note that data points at a given wavelength band were slightly shifted randomly 
in the horizontal direction for clarity. See the text for details. 


5 PROPERTIES OF THE CIRCUMSTELLAR DISKS 

The near- to mid-IR excesses above the stellar photosphere level probe the disk properties, such as the 
disk luminosity (basically an integral of the IR excesses) and disk geometry (e.g. Dullemond et al. 2007; 
Espaillat et al. 2013; Hughes et al. 2010; Kim et al. 2009; Merfn et al. 2010). In particular, disk flaring 
(e.g. Kenyon & Hartmann 1987) and radius of the inner disk edge are the two primary disk geometry 
parameters that shape the SED of IR excesses. As the disk evolves, dust grains in the inner circumstellar 
disks may gradually settle down (e.g. Dullemond & Dominik 2005) or be cleared out dynamically 
(Lubow & D’Angelo 2006) or through photoevaporation (Alexander et al. 2006a), which leads to a 
suppression of emission excesses progressively from near- to mid-IR wavelengths. Eeatures of the IR 
SEDs that are closely related to the disk clearing and flaring include the longest measured wavelength 
Atumoff shortward of which the emission is consistent with being purely from stellar photosphere, and 
the spectral index aexcess at A > Atumoff (e-g- Cieza et al. 2007; Harvey et al. 2007; Merin et al. 2008). 

5.1 aexcess VS. Atumoff 

Atumoff is closely related to the physical scales of the inward disk truncation or clearing radius (e.g. 
Calvet et al. 2002; Rice et al. 2003), and aexcess is related to both the inward disk clearing and disk 
flaring which in turn affect the disk temperature gradients. In particular, for an optically thick disk, a 
larger spectral index corresponds to a shallower temperature gradient (e.g. Beckwith et al. 1990). By 
comparing the observed SED of each YSO with the best-fit emergent stellar fluxes (be corrected for 
interstellar extinction), we determined the turnoff wavelength band Atumoff, longward of which > 3(7 
excesses above the stellar photosphere level were observed, and calculated aexcess for wavelength ranges 
longward of Atumoff- Previous studies of YSO IR spectral indices did not exclude the contribution of 
direct stellar photosphere emission. In this work, we focus on aexcess determined for the photosphere- 
subtracted IR SEDs in order to investigate the disk properties. In Table|2]we list spectral indices deter¬ 
mined for both the photosphere-included SEDs (atumoff) and photosphere-subtracted SEDs (aexcess) at 
A > Atumoff- The distribution on Atumoff vs. aexcess diagram for the subsamples of Stages I, II and III 
YSOs are shown separately in Eigure|9l 

The median aexcess at each Atumoff is also indicated as red triangles in Eigure|9] The majority of 
Stage I YSOs have aexcess ^ 0.0, whereas the majority of Stage HI YSOs have aexcess ^ 0.0. Compared 
to the Stages I and HI YSOs, the Stage H YSOs have a larger range of aexcess from ^ — 1 to 3. The 





14 


H.-X. Zhang, Y. Gao, & M. Fang et al. 



Stage I 



Stage II 

Stage III 

10“ 

A f i 4 

■ J 

1 

A 1 


n 

1 

10-* 

□ A 

r A S 1 


A 

1 





A 

A A 

10^“ 


i 



.. 

□< 

<]’ 

< 

_1 .. 


J H Ks IRl IR2 IR3 IR4 W3 

J H Kg 

IRl IR2 IR3 

IR4^ W3 

J H Kg IRl IR2 IRS IR4 ^3 


10 " 


10 ^ 


10 " 


10 * 


10 “ 


10 * 


'^turnoff 


'^turnoff (/^ni) 


'^turnoff (/^ni.) 


Fig. 10 Distributions of the fractional dust luminosity Ldust/f^* vs. Atumoff- YSOs of dif¬ 
ferent evolutionary stages are plotted separately in different panels. The median Tdust/f^* at 
different Atumoff is represented as red triangles. 


median ttexcess gradually increases with increasing Atumoff for both Stages I and II YSOs. No obvious 
trend of the median ctexcess with Atumoff is found for the Stage III YSOs. In addition, there is a hint 
that the standard deviation of ctexcess increases with increasing Atumoff for the Stage II YSOs which 
have the largest sample size. In particular, the standard deviations of ciexcess for the Stage II YSOs of 
different turnoff wavelengths at Atumoff £ IR2 and > IR3 are ^ 0.4 and 0.8 respectively. A smaller 
spread of atumoff at shorter Atumoff has been observed before (e.g. Cieza et al. 2007, Merln et al. 2008). 
Cieza et al. (2007) found that all the known Classical T Tauri stars (CTTs), which are defined by having 
relatively strong nebular emission lines and thus being actively accreting, cluster around atumoff ~ —1.0 
and Atumoff ^ Kg, whereas the Weak-line T Tauri stars (WTTs) exhibit a much larger spread of atumoff 
and Atumoff ■ 

5.2 Fractional Dust Luminosity vs. Atumoff 

The ratio of the circumstellar dust luminosity Ldust to stellar luminosity L*, which is also known as 
the fractional dust luminosity, was found to be correlated with the disk accretion activity (e.g. Kenyon 
& Hartmann 1995; Muzerolle et al. 2003). In particular, for mildly flared dusty disks, Ldust/i* ^ 0.1- 
0.2 cannot be simply explained by dust reprocessing of stellar radiation alone (Kenyon & Hartmann 
1995) but indicates that a significant amount of IR excesses may be contributed by self-radiation of an 
actively accreting disk, whereas YSOs with 0.001 ^ Ldust/7^* ^0.1 are expected to be mostly evolved 
objects with weaker or no observable accretion activity (e.g. Cieza et al. 2007). Moreover, most gas- 
poor debris disks (systems which are dominated by second-generation dust produced by the collision of 
planetesimals) were found to have Ldust/L* well below 0.001 (e.g. Currie & Kenyon 2009; Eiroa et al. 
2013; Matthews et al. 2014; Su et al. 2006; Trilling et al. 2008). 

We determined Ldust as an integral of the R06 model SED of the circumstellar dust (disk-tenvelope) 
that best fits the emergent IR excess emission, and L* as (i?*/i?Q)^(L^/T0)^, where i?* and are the 
stellar radius and effective temperature. The distribution of YSOs of different Stages on the Ldust/L* 
vs. Atumoff diagram is shown in EigurefTol The Stages I and HI YSOs are well separated at Ldust/L* 
0.1, whereas the Stage H YSOs have a range of Ldust/L* from ^ 0.01 to 1. Moreover, there is a general 
trend that the median Ldust/L* decreases with increasing Atumoff for YSOs of different evolutionary 
stages, pointing to an inside-out disk clearing of at least the small dust grains. 
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Fig. 11 Distribution of disk inner radius Rm vs. ctexcess for different Atumoff- R'm is normal¬ 
ized to the dust sublimation radius i?sub- YSOs with different ranges of Ldust/1^* are plotted 
with different symbols, as indicated in the top left panel. Distribution for the full sample is 
shown in the bottom right panel, where objects with 70/im detections are shown as red open 
circles. 


5.3 Disk Inner Radius vs. ctexcess 

The inner radius of the dusty disk determines the highest temperature of dust grains orbiting around 
the central stellar source (e.g. Backman & Paresce 1993), and thus can affect aexcess- Figure [TT] shows 
the relation between ctexcess and i?in/f?sub with different Atumoff. where i?in is the disk inner radius, 
and i?sub is the dust sublimation radius by assuming a sublimation temperature of 1600 K (R06). The 
bottom right panel of Figure [TT] shows the corresponding distribution for the full sample. 

There is a positive correlation between i?in/f?sub and aexcess at i?in/f?sub ^ 10 and aexcess ^ 0.0, 
irrespective of Ldust/1^* and Atumoff- A similar trend (not shown in the paper) also exists between i?in 
and aexcess at Rin > 0.5 AU and aexcess ^ 0.0. We note that a positive correlation was also found 
between disk inner radii (or hole radii) and disk masses for 35 c2d YSOs by Merin et al. (2010). 

5.4 Disk Flaring vs. aexcess 

Compared to a completely flat disk geometry, a flaring geometry increases the disk area that con¬ 
fronts the stellar radiation at large radii, and thus enhances the mid- to far-IR emission (e.g. Kenyon 
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Fig. 12 Distribution of disk flaring power /3 vs. ttexcess for different Atumoff- As in Figure 
m YSOs with different ranges of Ldust/i* are plotted with different symbols, as indicated 
in the top left panel. Distribution of the full sample is shown in the bottom right panel, where 
objects with 70pm detections are shown as red open circles. 
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& Hartmann 1987; Chiang & Goldreich 1997). The disk flaring power /3 describes the radial gradient 
of the disk scale height h, i.e. h{r) oc r^, where r is the cylindrical radius along the disk. Relationship 
between /3 and aexcess for our YSOs is shown in Figure [12] While no significant correlation between /3 
and aexcess was found for the overall sample, a majority of the disks with aexcess < 0.0 follow a trend 
that aexcess increases with /3, suggesting that the lack of a correlation between i?in/7?sub and aexcess for 
small disk inner radii can be in part attributed to the disk flaring. 

5.5 Discussion 

A variety of physical mechanisms have been invoked to explain the circumstellar disk evolution 
and clearing processes (e.g. Henning & Meeus 2011; Williams & Cieza 2011). The few commonly- 
considered mechanisms include viscous disk accretion (e.g. Hartmann et al. 1998; Lynden-Bell & 
Pringle 1974; Shakura & Sunyaev 1973), grain growth and dust settling (e.g. Dullemond & Dominik 
2005; Tanaka et al. 2005), photoevaporative dispersal (e.g. Alexander et al. 2006a,b; Gorti & Hollenbach 
2009; Hollenbach et al. 1994; Shu et al. 1993;) and dynamical clearing by companion stars or planets 
(e.g. Artymowicz & Lubow 1994; Kley & Nelson 2012; Lubow & D’Angelo 2006; Zhu et al. 2012). 
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While all of these proposed processes may operate simultaneously, it is important to probe the dominant 
process(es) at different disk evolution stages. 

5.5.1 From Oexcess to Disk Geometry 

The near- to mid-lR aexcess is primarily affected by the inner disk clearing and outer disk flaring. In par¬ 
ticular, the edge region of the optically thick inner disk, which is determined by either dust sublimation 
or some clearing processes, is illuminated directly by the stellar irradiation and thus contributes most 
of the hot dust emission excess, with the irradiation peak of this inner edge being shifted from near- to 
mid-lR as the disk is progressively cleared inside-out. In addition, as the disk evolves, dust settling or 
other clearing processes may result in a gradual reduction of disk flaring, which would in turn reduce 
the disk area that intercepts the stellar radiation and thus suppress the reprocessed cooler dust emission. 
Therefore, a progressively increasing disk inner edge is expected to increase ctexcess, whereas a smaller 
flaring power in the outer disk can result in a smaller aexcess- 

Our results suggest that variation of aexcess above ^ 0.0 primarily reflects the variation of disk 
clearing radii, whereas variation of aexcess below ~ 0.0 is largely related to a variation of the disk flaring 
power. Disk flaring is only important in shaping the near-to mid-IR SEDs when Rin ^ 10xi?sub (> 0.5 
AU for our sample). The lack of correlation between aexcess and disk flaring power at Rin ^ 10xi?sub 
implies that either the outer disk geometry does not vary synchronously with the inside-out disk clearing 
processes or spectral slopes at A ^ 24pm are not sensitive to the outer disk flaring. The small sample 
size of our disks (especially those with aexcess > 0.0) with detection at 70pm which is more sensitive to 
the outer disk flaring than shorter wavelengths (e.g. Sicilia-Aguilar et al. 2015) makes it hard to ascertain 
whether or not the outer disk flaring decreases or increases as the disk is cleared from the inside out. 
Recent studies of transitional disks in several nearby star-forming regions by Howard et al. (2013) and 
Keane et al. (2014) found that the continuum normalized [O l] 63.18pm line luminosities, which traces 
the cool, outer disks, are suppressed by a factor of ^ 2 on average with respect to the classical full disks, 
and this suppression was attributed to reduction of either the outer disk flaring or gas-to-dust ratio. 

5.5.2 Probing Disk Dispersal Processes with Transitional Disks 

There may be a variety of evolutionary paths from the optically thick full disks to optically thin to debris 
disks. Distinguishing different disk dispersal processes is crucial for understanding how the planetary 
systems formed from protoplanetary disks. The partially-cleared transitional disks, which have little or 
no excess emission in the near-IR (^ 5pm) and thus optically thin inner opacity holes but a significant 
excess at longer wavelengths (e.g. Brown et al. 2007; Calvet et al. 2005; Strom et al. 1989; Skrutskie et 
al. 1990), provide a unique opportunity to probe different disk clearing mechanisms because different 
mechanisms are expected to result in very different IR spectral slopes, disk luminosities, and accretion 
activities in the short transitional stages (e.g. Alexander et al. 2014; Cieza et al. 2010; Najita, Strom & 
Muzerolle 2007). 

To open an inner opacity hole through photoevaporation, the disk viscous accretion rate has to fall 
below the photoevaporation rate (e.g. Alexander, Clarke & Pringle 2006a; Owen et al. 2010), and once 
this happens, the full disks of gas and small dust grains can be quickly dissipated from the inside out 
in ^ 0.1 Myr which is an order of magnitude shorter than the typical disk lifetime. Besides a low frac¬ 
tional disk luminosity and steep IR spectral slope (e.g. aexcess < 0.0), another important consequence 
from photoevaporative clearing is that little or no accretion is expected once an inner hole is opened. In 
contrast, dynamical clearing by giant planets may sustain a small but still considerable amount of disk 
accretion across the inner opacity hole and a relatively high outer disk masses and luminosities and thus 
rising mid- to far-IR SEDs (e.g. Alexander 2008; Najita et al. 2007). In contrary to both photoevapora¬ 
tion and dynamical clearing, the pure grain growth and dust settling processes can result in an efficient 
depletion of small grains (and thus suppression of near- to mid-IR emission) from the inside out over 
time scales much smaller than 0.1 Myr (e.g. Dullemond & Dominik 2005), with little direct influence 
on the accretion activity. 
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Fig. 13 ^excess IS plotted agaiiist the fractional disk luminosities Ldust/1^*- Disks with 
Aturnoff < and > IR3 (transitional disks) are plotted as red open circles and black filled 
squares respectively. The horizontal dashed line separates the sample into disks with 
Ldust/Li, > and <0.1, and the vertical dashed line separates the sample into disks with 
ctexcess > and < 0.0. Most accreting disks were found to have Ldust/T* >0.1. Transitional 
disks in the lower left part may be primarily cleared by photoevaporation, while those in the 
upper right part may be dynamically cleared by giant planets. 


All of our YSOs have Ldust/T* > 10“°, and 49 (23%) have Aturnoff > IR3 and thus can be 
classified as transitional disks. Recall that for our sample disks with Aturnoff > IR3 exhibit a remarkably 
higher median and larger scatter of ctexcess than those with Aturnoff < IR3 (Figure [O]). The fraction of 
transitional disks in our sample is slightly higher yet still comparable to previous studies of nearby star 
clusters or star-forming regions (e.g. Currie et al. 2009; Dahm & Carpenter 2009; Fang et al. 2009; 
Hernandez et al. 2007; Kim et al. 2009; Lada et al. 2006). The distribution of our sample on the aexcess 
vs. Ldust/T* plane is shown in Figure [13] where the transitional disks are plotted as black squares 
(filled for those with < 4000 K, and open for those with T* > 4000 K). Note that previous studies 
did not subtract the stellar photosphere emission for calculating the excess spectral index, which tends 
to underestimate the “genuine” ctexcess- 

As is shown in Figure [T3 the majority of the disks with Aturnoff < IR3 clustered toward the upper 
left corner, with Ldust/T* ^ 10“^ and aexcess ^ 0.0, whereas the disks with Aturnoff > IR3 seem to 
follow a sequence from the upper right to the lower left, with none of them having Ldust/T* >10“^ 
and aexcess < 0.0. Most of the objects around the upper left corner are expected to have accreting full 
disks, and they are clearly separated from the population of transitional disks in Figure [T^ A similar 
separation of transitional disks and full disks was also recently found by Sicilia-Aguilar et al. (2015) 
based on the relation between spectral indices and accretion rates. Among the objects with Aturnoff < 
IR3, 14 (7%) have aexcess < 0.0 and Ldust/T* < 0.003. These 7% objects are consistent with being the 
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so-called “anemic” (e.g. Lada et al. 2006) or “homologously depleted” (e.g. Currie & Sicilia-Aguilar 
2011) disks, which have detectable excess emission that decreases steadily at all wavelengths. 

Transitional disks toward the lower left corner of Figure[T3]may be more evolved than those toward 
the upper right. Among the 49 transitional disks, 41 (84%) have aexcess > 0.0 and 8 (16%) have 
ctexcess < 0.0. Observations of UV continuum or recombination emission lines for all of our sample 
will be necessary for obtaining ongoing disk accretion rate. The accretion activities are known to be 
closely connected to the disk global properties, such as disk luminosities, masses and dust settling. If 
we instead use Ldust/^* to approximately discriminate disks with and without accretion activity at 
a dividing value = 0.1, 17 (35%) of the 49 transitional disks have aexcess > 0.0 and Ldust/i* > 
0.1, which may indicate the possibility of dynamical clearing by giant planets; Among the 32 (65%) 
disks with Ldust/^* < 0.1, 8 have aexcess < 0.0 and 24 have aexcess > 0.0, the low Ldust/^* 
probably indicates that these 32 disks are primarily cleared by photoevaporation. None of our transitional 
disks have aexcess < 0.0 and Ldust/7j* > 0.1, so grain growth and dust settling alone are probably 
not important hole-opening mechanisms (Ceiza et al. 2010). Furthermore, our finding that the median 
aexcess of Stages I and II YSOs tend to increases with Atumoff also suggests that disk clearing is not 
primarily driven by grain growth which would otherwise results in a negative correlation between aexcess 
and Aturnoff (e.g. Dullemond & Dominik 2005). 

6 SUMMARY 

We have statistically explored the properties of the central stellar sources, the evolutionary stages, and 
the circumstellar disks for a sample of 211 Perseus YSOs by modeling the optical to mid-IR broadband 
SEDs with the R06 YSO evolution models. The median central stellar mass and age for the Perseus 
YSOs are ~ 0.3 Mq and ~ 3.1 Myr respectively based on the Siess et al. (2000) PMS evolutionary 
models. About 81% of our sample are classified as Stage II objects which are characterized by having 
optically thick disks, ^ 14% are classified as Stage I objects which are characterized by having signif¬ 
icant infalling envelopes, and the remaining 5% are classified as Stage III objects with optically thin 
disks. Our primary results are summarized as follows. 

- The evolutionary Stages as determined from the SED modeling have a general correspondence with 
the traditional spectral-indices-based Classes. In particular, ^ 90% of the Class II YSOs fall into 
the Stage II phase which is characterized by optically thick disks, and 75% of the Class I YSOs 
fall into the Stage I phase which is characterized by significant infalling envelopes. Nevertheless, 
relating the Classes III and Elat YSOs to specific evolutionary stages is uncertain. In particular, 
half of the Class III YSOs fall into the Stage II and the other half fall into the optically thin Stage 
III phase, and half of the Class Elat YSOs fall into the Stage I and the other half fall into the Stage 
II phase. 

- We determined the turnoff wave band (Atumoff) longward of which significant IR excesses with 
respect to the stellar photosphere level start to be observed and the excess spectral indices aexcess 
at A > Aturnoff- The median and standard deviation of aexcess of the Stages I and II YSOs tend to 
increase with Atumoff, especially at Atumoff > IRAC 5.8/rm. There is a general trend that the median 
fractional dust luminosity Ldust/E* decrease with increasing Atumoff, pointing to an inside-out disk 
clearing of small dust grains. We found a positive correlation between aexcess and disk inner radius 
Rin, and a lack of correlation between aexcess and disk flaring at aexcess ^0.0 and i?in ^ 10 x i?sub, 
which indicates that, firstly, the near- to mid-IR spectral slopes primarily reflect the progressive disk 
clearing from the inside out once Rin ^ 10xi?sub, secondly, the outer disk flaring either does not 
vary synchronously with the inner disk clearing processes or has little appreciable influence on the 
spectral slopes at wavelengths ^ 24^m. 

- About 23% (49) of our YSOs are classified as transitional disks, which have Atumoff > IRAC 
5.8^m and Ldust/E* > 10“^. By using the Ldust/E* to approximately discriminate disks with and 
without accretion activity at a dividing value of 0.1, 35% of the transitional disks have aexcess > 0.0 
and Ldust/T* > 0.1, implying the possibility of dynamical clearing by giant planets; 65% have 
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Ld-ast/Li, < 0.1, which is consistent with the expectation of photoevaporative clearing; None of 
our disks has ctexcess < 0.0 and idust/^* > 0.1, so grain growth and dust settling are probably 
not the driven mechanisms in disk clearing, in line with the trend that the median ttexcess increases, 
rather than decreases, with Atumoff ■ 

An indispensable diagnostic of the evolutionary stages of YSOs and their circumstellar disks is 
the current accretion rate, which is usually determined either from recombination lines or ultraviolet 
continuum excesses. Different disk clearing processes can lead to different disk accretion properties, 
the effect of which is especially prominent in transitional stages. Moreover, similar to many previous 
studies, our current work is heavily biased against the Stage III YSOs with optically thin or anemic 
disks. To understand the disk evolution and dispersal processes, a systematic census of Stage III YSOs 
and their disk accretion activity is imperative. Therefore, our future direction will include 1) a system¬ 
atic spectroscopic followup of our YSOs with the Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST; Cui et al. 2012) to place stringent constraints on the on-going accretion activity; 
2) a wide-held time-series optical photometry across the whole Perseus region for an unbiased census 
of Stage III disks with the PMO Xuyi 1.2-m Schmidt Telescope, in order to further probe the dominant 
disk dispersal mechanisms. 


Table 1: Optical Photometry of Perseus YSOs 


IDc2d 

R.A.(J2000) 

Decl.(J2000) 

9 

o's 

r 

0" 'p 

i 

O'* 


(deg) 

(deg) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

J band turnoff^ 

c2dJ032852.2+304506 

52.21736 

30.75154 

13.96 

0.06 

12.9 

0.06 



c2dJ032854.6+311651 

52.22763 

31.28086 



18.73 

0.08 

16.7 

0.06 

c2dJ032917.7+312245 

52.32366 

31.37917 

15.55 

0.06 

13.82 

0.05 



c2dJ033035.5+311559 

52.64781 

31.26627 

18.78 

0.07 

18.0 

0.07 

17.05 

0.06 

c2dJ033330.4+311051 

53.37669 

31.18071 

15.99 

0.05 

13.98 

0.06 



c2dJ034157.4+314837 

55.48934 

31.81021 

18.38 

0.06 

16.22 

0.06 

14.7 

0.06 

c2dJ034344.5+314309 

55.93535 

31.71926 

18.69 

0.06 

17.19 

0.06 

14.52 

0.06 

c2dJ034413.0+320135 

56.05407 

32.02652 





18.99 

0.13 

c2dJ034441.2+321010 

56.17156 

32.16944 





17.4 

0.07 

H band turnoff 

c2dJ032519.5+303424 

51.33134 

30.57338 

19.24 

0.07 

18.75 

0.09 

18.34 

0.06 

c2dJ033037.0+303128 

52.65402 

30.52437 

20.41 

0.11 

17.82 

0.06 

15.82 

0.06 

c2dJ033044.0+303247 

52.68326 

30.54639 

16.13 

0.07 

14.55 

0.07 

13.21 

0.08 

c2dJ033118.3+304940 

52.82625 

30.82765 

18.11 

0.06 

15.31 

0.06 

14.48 

0.06 

c2dJ033312.8+312124 

53.30349 

31.35673 



19.61 

0.13 

17.88 

0.06 

c2dJ033341.3+311341 

53.42204 

31.22806 

20.9 

0.32 

18.52 

0.07 

17.73 

0.06 

c2dJ034109.1+314438 

55.28804 

31.74386 

17.63 

0.06 

15.18 

0.06 

14.28 

0.06 

c2dJ034255.9+315842 

55.73312 

31.97834 

14.69 

0.06 

13.45 

0.06 



c2dJ034426.7+320820 

56.11124 

32.13898 

19.45 

0.12 

18.34 

0.08 

16.7 

0.06 

c2dJ034431.1+321848 

56.12973 

32.31347 



19.86 

0.37 

18.09 

0.08 

c2dJ034437.9+320804 

56.15785 

32.13448 

18.1 

0.06 

16.13 

0.06 

14.87 

0.06 

c2dJ034516.3+320620 

56.31809 

32.10559 

18.33 

0.06 

16.19 

0.06 

14.85 

0.06 

c2dJ034520.5+320634 

56.33525 

32.10958 

17.57 

0.06 

15.39 

0.06 

14.1 

0.06 

Ks band turnoff 

c2dJ032741.5+302017 

51.92281 

30.33799 

16.46 

0.06 

14.74 

0.06 



c2dJ032800.1+300847 

52.00038 

30.1464 

19.82 

0.08 

17.32 

0.06 

15.13 

0.07 

c2dJ032847.6+312406 

52.19853 

31.40168 





19.64 

0.11 

c2dJ032850.6+304245 

52.2109 

30.7124 





20.31 

0.28 

c2dJ032851.0+311818 

52.21262 

31.30513 

18.49 

0.07 

15.81 

0.07 

14.0 

0.07 

c2dJ032851.2+311955 

52.21335 

31.3319 

18.19 

0.06 

15.94 

0.06 

14.7 

0.06 

c2dJ032859.6+312147 

52.24817 

31.36296 

18.42 

0.06 

16.13 

0.06 

15.59 

0.06 

c2dJ032903.8+311604 

52.26574 

31.26773 

20.12 

0.13 

17.61 

0.07 

16.01 

0.06 

c2dJ032903.9+305630 

52.26613 

30.9416 





18.87 

0.07 

c2dJ032903.9+312149 

52.26614 

31.3635 

17.15 

0.07 

15.28 

0.08 

14.18 

0.07 

c2dJ032909.0+312624 

52.28738 

31.43997 





19.95 

0.17 


The YSOs start exhibiting significant (3c7) IR excesses above the photosphere level longwai'd of the turnoff wavebands. 
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IDc2d 

R.A.(J2000) 

Decl.(J2000) 

9 


r 

(7 r 

i 

o-i 


(deg) 

(deg) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

c2dJ032910.8+311643 

52.29515 

31.27849 



19.79 

0.16 

18.71 

0.07 

c2dJ032913.1+312253 

52.30474 

31.38134 



19.71 

0.19 

17.3 

0.07 

c2dJ032921.9+311536 

52.34115 

31.26005 

15.29 

0.05 

14.96 

0.06 

13.49 

0.06 

c2dJ032923.2+312030 

52.34653 

31.34173 

18.59 

0.06 

16.9 

0.06 

15.14 

0.06 

c2dJ032932.6+312437 

52.38573 

31.41025 



17.79 

0.06 

16.65 

0.06 

c2dJ033001.9+303529 

52.5078 

30.59145 



19.09 

0.1 

18.21 

0.08 

c2dJ033035.9+303024 

52.64968 

30.50678 



12.07 

0.07 



c2dJ033038.2+303212 

52.65919 

30.53665 



19.7 

0.2 

18.44 

0.07 

c2dJ033052.5+305418 

52.71878 

30.90494 



19.05 

0.12 

17.14 

0.06 

c2dJ033114.7+304955 

52.81127 

30.83206 

18.53 

0.06 

16.91 

0.06 

15.8 

0.06 

c2dJ033142.4+310625 

52.92668 

31.10691 

19.8 

0.1 

17.66 

0.06 

15.8 

0.06 

c2dJ033233.0+310222 

53.13745 

31.03935 

20.13 

0.11 

17.84 

0.06 

15.38 

0.06 

c2dJ033234.0+310056 

53.14185 

31.01549 

18.02 

0.06 

15.79 

0.06 

14.34 

0.06 

c2dJ033241.7+311046 

53.17377 

31.17953 

20.63 

0.2 

18.35 

0.07 

16.65 

0.06 

c2dJ033401.7+311440 

53.50692 

31.24438 

17.86 

0.06 

15.65 

0.06 

13.95 

0.06 

c2dJ033915.8+312431 

54.81587 

31.40854 





18.77 

0.07 

c2dJ034119.2+320204 

55.32994 

32.03438 

20.07 

0.13 

17.76 

0.07 

16.15 

0.06 

c2dJ034155.7+314811 

55.48214 

31.80318 





20.11 

0.12 

c2dJ034157.8+314801 

55.49064 

31.80023 



19.16 

0.08 

16.9 

0.06 

c2dJ034219.3+314327 

55.5803 

31.72415 



19.56 

0.1 

17.44 

0.06 

c2dJ034232.9+314221 

55.63711 

31.70572 



17.63 

0.06 

16.16 

0.06 

c2dJ034322.2+314614 

55.84257 

31.77045 



19.32 

0.11 

17.45 

0.06 

c2dJ034328.2+320159 

55.86753 

32.0331 

16.88 

0.06 

15.68 

0.06 

14.52 

0.06 

c2dJ034355.2+315532 

55.98018 

31.92559 





18.9 

0.08 

c2dJ034356.0+320213 

55.98346 

32.03702 

20.51 

0.17 

18.09 

0.07 

16.35 

0.06 

c2dJ034358.6+321728 

55.99406 

32.29097 



16.32 

0.1 

14.93 

0.07 

c2dJ034358.9+321127 

55.99549 

32.19088 

18.06 

0.06 

16.36 

0.06 

15.3 

0.05 

c2dJ034359.9+320441 

55.9995 

32.07817 





17.67 

0.06 

c2dJ034406.0+321532 

56.02504 

32.25892 





17.74 

0.07 

c2dJ034406.8+320754 

56.02833 

32.13167 

19.46 

0.1 

17.81 

0.07 

16.1 

0.06 

c2dJ034407.5+320409 

56.03132 

32.0691 

20.37 

0.17 

18.28 

0.08 

16.48 

0.06 

c2dJ034411.6+320313 

56.04844 

32.05364 



19.03 

0.17 

17.22 

0.06 

c2dJ034418.6+321253 

56.07747 

32.21475 



18.66 

0.1 

19.08 

0.13 

c2dJ034421.6+321038 

56.0901 

32.17713 

18.99 

0.08 

16.85 

0.06 

15.72 

0.06 

c2dJ034422.3+321201 

56.09307 

32.20019 

18.31 

0.06 

16.67 

0.06 

15.26 

0.06 

c2dJ034425.5+321131 

56.10633 

32.192 

19.28 

0.08 

17.22 

0.06 

15.63 

0.06 

c2dJ034427.3+321421 

56.11359 

32.23915 

19.42 

0.08 



16.0 

0.06 

c2dJ034431.4+320014 

56.13069 

32.00394 





18.56 

0.08 

c2dJ034435.7+320304 

56.1487 

32.05097 

20.47 

0.2 

18.88 

0.13 

17.55 

0.06 

c2dJ034438.5+320736 

56.16024 

32.12659 

16.84 

0.06 

15.25 

0.06 

14.38 

0.06 

c2dJ034438.5+320801 

56.1606 

32.13351 

18.92 

0.07 

16.9 

0.06 

15.21 

0.06 

c2dJ034444.7+320402 

56.18633 

32.06736 

17.54 

0.06 

15.42 

0.06 

14.06 

0.06 

c2dJ034452.0+322625 

56.21668 

32.4404 



18.83 

0.1 

16.92 

0.06 

c2dJ034452.1+315825 

56.21689 

31.97367 





17.71 

0.06 

c2dJ034525.1+320930 

56.35479 

32.15842 

19.0 

0.07 

16.85 

0.06 

15.06 

0.06 

c2dJ034536.8+322557 

56.40347 

32.43251 

15.96 

0.06 

13.79 

0.05 



c2dJ034548.3+322412 

56.45111 

32.40334 





10.71 

0.07 

c2dJ034558.2+322647 

56.49269 

32.44653 

20.7 

0.29 

18.72 

0.1 

16.73 

0.06 

IRl band turnoff 

c2dJ032747.7+301205 

51.94864 

30.20126 



18.17 

0.07 

16.14 

0.07 

c2dJ032834.5+310051 

52.1437 

31.01419 





18.38 

0.11 

c2dJ032842.4+302953 

52.17673 

30.4981 

17.97 

0.06 

16.28 

0.06 

14.53 

0.06 

c2dJ032844.1+312053 

52.18372 

31.34799 





18.34 

0.06 

c2dJ032846.2+311638 

52.19252 

31.27734 

16.63 

0.06 

14.93 

0.06 

13.48 

0.06 

c2dJ032847.8+311655 

52.19933 

31.28196 

20.02 

0.12 

19.03 

0.08 

16.6 

0.06 

c2dJ032852.2+312245 

52.2174 

31.37924 

17.9 

0.06 

15.99 

0.06 

14.64 

0.06 

c2dJ032856.6+311836 

52.23602 

31.30987 



17.65 

0.07 

16.28 

0.06 

c2dJ032857.0+311622 

52.23736 

31.27285 





18.12 

0.07 

c2dJ032903.1+312238 

52.26311 

31.37723 

20.11 

0.14 

18.47 

0.07 

17.06 

0.06 

c2dJ032904.1+305613 

52.26716 

30.9369 

20.58 

0.16 

19.41 

0.12 

17.06 

0.06 

c2dJ032918.7+312325 

52.32808 

31.39038 

16.0 

0.05 

14.15 

0.06 

13.57 

0.06 

c2dJ032920.4+311834 

52.33515 

31.3095 





18.72 

0.13 
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IDc2d 

R.A.(J2000) 

Decl.(J2000) 

9 


r 

(7 r 

i 

o-i 


(deg) 

(deg) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

c2dJ032930.4+311903 

52.37668 

31.31759 

17.57 

0.06 

16.05 

0.06 

14.69 

0.06 

c2dJ032932.9+312713 

52.387 

31.45349 

19.35 

0.07 

17.67 

0.06 

15.93 

0.06 

c2dJ032937.7+312202 

52.40723 

31.36735 



19.71 

0.12 

17.3 

0.06 

c2dJ032954.0+312053 

52.47518 

31.34803 

17.81 

0.06 

15.8 

0.06 

14.72 

0.06 

c2dJ033024.1+311404 

52.60034 

31.23454 

19.3 

0.07 

18.0 

0.06 

16.18 

0.06 

c2dJ033110.7+304941 

52.79451 

30.82795 

18.55 

0.06 

15.52 

0.06 

15.05 

0.06 

c2dJ033430.8+311324 

53.62826 

31.22343 



18.45 

0.07 

16.38 

0.06 

c2dJ033449.8+311550 

53.70768 

31.26396 

15.85 

0.06 

14.11 

0.06 



c2dJ034001.5+311017 

55.00621 

31.17147 



19.56 

0.15 

17.65 

0.06 

c2dJ034201.0+314913 

55.50422 

31.82038 





19.16 

0.07 

c2dJ034204.3+314712 

55.51807 

31.78655 





18.05 

0.06 

c2dJ034220.3+320531 

55.58467 

32.09195 

20.03 

0.12 

17.57 

0.06 

16.14 

0.06 

c2dJ034232.1+315250 

55.63377 

31.88043 





19.92 

0.14 

c2dJ034249.2+315011 

55.70492 

31.83643 



19.47 

0.1 

17.11 

0.06 

c2dJ034313.7+320045 

55.80708 

32.01254 





18.71 

0.08 

c2dJ034323.6+321226 

55.84821 

32.20718 



19.32 

0.12 

17.23 

0.06 

c2dJ034329.4+315219 

55.87265 

31.87207 





19.38 

0.1 

c2dJ034345.2+320359 

55.9382 

32.06628 





18.51 

0.1 

c2dJ034348.8+321552 

55.95345 

32.26431 

19.57 

0.13 

17.63 

0.07 

15.95 

0.06 

c2dJ034355.3+320753 

55.98033 

32.13147 



19.01 

0.13 

17.04 

0.06 

c2dJ034359.1+321421 

55.99624 

32.23923 

20.15 

0.13 

17.61 

0.07 

16.73 

0.06 

c2dJ034401.6+322359 

56.00656 

32.39968 



18.68 

0.11 

16.73 

0.06 

c2dJ034402.9+315228 

56.01215 

31.87437 





18.9 

0.09 

c2dJ034418.2+320457 

56.0757 

32.08249 



18.09 

0.07 

16.4 

0.06 

c2dJ034425.5+320617 

56.10645 

32.10476 



18.47 

0.09 

16.77 

0.06 

c2dJ034426.0+320430 

56.10848 

32.07512 

16.19 

0.06 

14.19 

0.06 

13.33 

0.06 

c2dJ034427.9+322719 

56.11625 

32.45525 

19.42 

0.08 

17.99 

0.07 

16.06 

0.06 

c2dJ034428.5+315954 

56.1188 

31.99833 

18.65 

0.06 

16.69 

0.06 

15.03 

0.06 

c2dJ034429.8+320055 

56.12418 

32.01516 





17.49 

0.06 

c2dJ034432.0+321144 

56.1335 

32.19548 

15.47 

0.06 

13.8 

0.06 

12.77 

0.06 

c2dJ034433.8+315830 

56.1408 

31.97506 



18.75 

0.12 

16.88 

0.06 

c2dJ034435.0+321531 

56.1458 

32.25865 



18.36 

0.08 

16.5 

0.06 

c2dJ034435.5+320856 

56.14779 

32.14897 





16.9 

0.06 

c2dJ034437.0+320645 

56.15399 

32.11256 



12.16 

0.06 



c2dJ034437.4+320901 

56.1559 

32.15024 



16.75 

0.09 

15.54 

0.06 

c2dJ034438.0+320330 

56.15825 

32.05825 

17.72 

0.06 

15.43 

0.06 

14.14 

0.06 

c2dJ034439.8+321804 

56.16583 

32.30112 

19.6 

0.09 

17.12 

0.06 

15.44 

0.06 

c2dJ034440.2+320933 

56.16771 

32.15917 





17.01 

0.08 

c2dJ034442.6+321002 

56.17741 

32.16735 





17.97 

0.07 

c2dJ034443.1+313734 

56.17942 

31.62603 





19.06 

0.08 

c2dJ034443.8+321030 

56.18241 

32.1751 

20.86 

0.24 

18.55 

0.1 

15.89 

0.06 

c2dJ034450.4+315236 

56.20979 

31.87667 





18.27 

0.09 

c2dJ034456.1+320915 

56.23394 

32.15422 

17.08 

0.06 

14.75 

0.06 

13.77 

0.05 

c2dJ034517.8+321206 

56.32426 

32.20162 



18.78 

0.09 

16.78 

0.06 

c2dJ034529.7+315920 

56.37382 

31.98881 





18.53 

0.08 

c2dJ034533.5+314555 

56.38945 

31.76536 





19.73 

0.13 

c2dJ034535.6+315954 

56.39849 

31.99845 





17.88 

0.07 

c2dJ034657.4+324917 

56.7391 

32.8215 



18.7 

0.09 

16.86 

0.06 

1R2 band turnoff 

c2dJ032851.1+311632 

52.21281 

31.27566 

19.99 

0.08 

18.22 

0.07 

16.27 

0.06 

c2dJ032852.2+311547 

52.2173 

31.26307 



19.19 

0.1 

16.83 

0.06 

c2dJ032852.9+311626 

52.22052 

31.274 

20.38 

0.15 

18.47 

0.07 

16.53 

0.06 

c2dJ032909.5+312721 

52.28954 

31.45581 



20.12 

0.23 

17.96 

0.06 

c2dJ032917.8+311948 

52.32406 

31.33001 





19.35 

0.08 

c2dJ032921.6+312110 

52.33988 

31.35287 

18.91 

0.07 

16.85 

0.06 

15.24 

0.06 

c2dJ032923.2+312653 

52.34687 

31.44808 

20.04 

0.11 

18.57 

0.07 

16.69 

0.06 

c2dJ032928.9+305842 

52.37034 

30.9783 

18.75 

0.06 

17.2 

0.06 

15.42 

0.06 

c2dJ032929.8+312103 

52.37415 

31.35072 

19.65 

0.1 

17.72 

0.06 

15.83 

0.06 

c2dJ032937.6+310249 

52.40678 

31.04699 





17.85 

0.06 

c2dJ033120.1+304918 

52.83379 

30.82157 



19.16 

0.09 

17.07 

0.06 

c2dJ033346.9+305350 

53.44552 

30.89726 

20.43 

0.17 

18.4 

0.07 

16.55 

0.06 

c2dJ034336.0+315009 

55.90007 

31.83583 





19.66 

0.09 
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IDc2d 

R.A.(J2000) 

Decl.(J2000) 

9 


r 

(7 r 

i 

o-i 


(deg) 

(deg) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

c2dJ034346.5+321106 

55.94357 

32.18498 





18.76 

0.09 

c2dJ034347.6+320903 

55.94853 

32.1507 





18.54 

0.09 

c2dJ034415.8+315937 

56.06598 

31.99354 





17.52 

0.06 

c2dJ034421.3+321156 

56.08878 

32.19897 

19.09 

0.07 

16.94 

0.06 

15.49 

0.06 

c2dJ034427.2+322029 

56.11342 

32.34133 





17.01 

0.06 

c2dJ034430.8+320956 

56.12848 

32.16547 

11.93 

0.08 



11.02 

0.15 

c2dJ034658.5+324659 

56.74379 

32.78303 



19.22 

0.13 

17.05 

0.06 

IR3 band turnoff 

c2dJ032858.1+311804 

52.24213 

31.30102 

19.06 

0.07 

16.92 

0.06 

15.65 

0.06 

c2dJ032908.0+312251 

52.28315 

31.38095 



19.26 

0.11 

17.22 

0.08 

c2dJ032912.9+312329 

52.30381 

31.39147 



18.36 

0.09 

16.65 

0.06 

c2dJ032916.8+312325 

52.32013 

31.39031 





18.22 

0.06 

c2dJ032926.8+312648 

52.36172 

31.44654 

15.81 

0.05 

14.11 

0.06 



c2dJ032929.3+311835 

52.37198 

31.30963 

19.27 

0.07 

17.01 

0.06 

15.6 

0.06 

c2dJ034233.1+315215 

55.63803 

31.87075 





18.86 

0.08 

c2dJ034234.2+315101 

55.64244 

31.85028 





18.23 

0.07 

c2dJ034250.9+314045 

55.71208 

31.67921 





18.44 

0.06 

c2dJ034301.9+314436 

55.75807 

31.74322 





19.4 

0.08 

c2dJ034308.7+315139 

55.78628 

31.86072 





19.06 

0.08 

c2dJ034344.6+320818 

55.93594 

32.13827 

17.3 

0.06 

15.46 

0.06 

14.46 

0.05 

c2dJ034410.1+320405 

56.0422 

32.06792 



19.2 

0.16 

16.93 

0.06 

c2dJ034415.2+321942 

56.06348 

32.32838 





17.6 

0.08 

c2dJ034418.2+320959 

56.07588 

32.16648 



19.23 

0.16 

17.14 

0.06 

c2dJ034422.3+320543 

56.09287 

32.09521 

18.94 

0.07 

16.87 

0.06 

15.61 

0.06 

c2dJ034422.6+320154 

56.09409 

32.03157 

18.87 

0.07 

16.75 

0.06 

15.05 

0.06 

c2dJ034425.7+321549 

56.10713 

32.26367 





18.83 

0.14 

c2dJ034429.2+320116 

56.1218 

32.02103 





18.9 

0.1 

c2dJ034429.7+321040 

56.12391 

32.17772 

17.39 

0.11 

15.87 

0.07 

14.83 

0.06 

c2dJ034434.1+321636 

56.14225 

32.2766 





17.57 

0.06 

c2dJ034434.8+315655 

56.14503 

31.94866 

19.52 

0.09 

17.47 

0.07 

15.83 

0.06 

c2dJ034437.4+321224 

56.15584 

32.20671 



18.47 

0.09 

16.92 

0.06 

c2dJ034438.0+321137 

56.15838 

32.19361 



18.3 

0.08 

16.45 

0.06 

c2dJ034439.0+320320 

56.16238 

32.05547 



19.05 

0.18 

17.39 

0.06 

c2dJ034439.2+322009 

56.16331 

32.3358 

19.89 

0.11 

17.87 

0.06 

15.83 

0.06 

c2dJ034441.7+321202 

56.17392 

32.20062 





16.98 

0.06 

c2dJ034442.6+320619 

56.17758 

32.10541 

17.01 

0.06 

15.4 

0.06 

14.58 

0.06 

c2dJ034443.0+321560 

56.17929 

32.26656 





18.15 

0.07 

c2dJ034444.6+320813 

56.18579 

32.13681 

19.53 

0.13 

17.8 

0.06 

16.16 

0.06 

c2dJ034457.9+320402 

56.24106 

32.0671 



18.63 

0.1 

16.56 

0.06 

c2dJ034460.0+322233 

56.24997 

32.37576 





17.76 

0.06 

c2dJ034501.4+320502 

56.25595 

32.08382 

13.76 

0.06 

12.79 

0.05 



c2dJ034504.7+321501 

56.2694 

32.2503 



19.27 

0.15 

17.14 

0.06 

c2dJ034513.5+322435 

56.30627 

32.40966 





17.89 

0.06 

IR4 band turnoff 

c2dJ032854.1+311654 

52.22537 

31.28172 



18.89 

0.08 

16.59 

0.06 

c2dJ033027.1+302830 

52.61309 

30.47493 





18.39 

0.1 

c2dJ034227.1+314433 

55.613 

31.74247 



19.17 

0.11 

17.33 

0.06 

c2dJ034254.7+314345 

55.72778 

31.72924 

16.22 

0.05 

14.51 

0.06 

13.59 

0.06 

c2dJ034306.8+314821 

55.77822 

31.80569 

20.05 

0.12 

17.63 

0.06 

15.63 

0.06 

c2dJ034419.1+320931 

56.07973 

32.15869 





11.3 

0.08 

c2dJ034421.6+321510 

56.08987 

32.25271 

18.57 

0.1 

17.62 

0.07 

15.83 

0.06 

c2dJ034431.5+320845 

56.13145 

32.14581 

14.71 

0.06 

13.4 

0.06 

12.76 

0.06 

c2dJ034456.8+315411 

56.23684 

31.90317 





17.45 

0.06 

c2dJ034507.6+321028 

56.28182 

32.17441 



12.11 

0.05 



W3 band turnoff 

c2dJ032916.7+311618 

52.31955 

31.27171 

17.17 

0.06 

15.52 

0.06 

13.97 

0.06 

c2dJ033026.0+310218 

52.60821 

31.03831 

15.27 

0.05 

13.72 

0.05 



c2dJ033351.1+311228 

53.46281 

31.20772 

19.35 

0.07 

17.03 

0.06 

14.77 

0.06 

c2dJ034011.8+315523 

55.04929 

31.92315 

19.42 

0.07 

17.39 

0.06 

16.17 

0.06 


Notes: All the photometry was calibrated against the SDSS DR8 catalog. 
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Table 2: SED Fitting Results of Perseus YSOs 


IDc2d 

lOjnodel 

xi Av 

D 

M* 

Age* L* 

-^dust 

-Rin/-f^sub 

Stage QturnofF (^excess 

( 1 ) 

( 2 ) 

(mag) 

(3) (4) 

(kpc) 

(5) 

(AIq) 

( 6 ) 

( 10 ^ yr) (IO-^Lq) 

(7) ( 8 ) 

(9) 

( 10 ) 

(11) (12) (13) 

J band turnoff^ 


c2dJ032852.2+304506 

3000104 

2.0 

3.76 

0.28 

2 36'^^ "^^ 

Z..JU_Q 22 

76.39:^-;- 

5439.09 

173.77 

^ Y+u. / 

^•^-0.7 

II 

-0.75 

-0.63 

c2dJ032854.6+311651 

3013176 

0.4 

3.04 

0.23 

A 14+0.23 
^•^^-0.04 

25.55127-4 

8.72 

1.13 

^•^-3.2 

II 

-0.57 

-0.54 

c2dJ032917.7+312245 

3005510 

0.5 

5.7 

0.26 

'y 'yy + 0.93 

/o.^z_gg 3 

2637.26 

263.42 

1 1+0.5 
^•^- 0.1 

II 

-0.95 

- 0.8 

c2dJ033035.5+311559 

3018343 

6.7 

0.79 

0.24 

1 02 ^^®^ 

Q Qn+31.75 
®-^^_8.84 

2475.96 

1910.17 

83 0+237.7 
o-^-^_82.9 

I 

-0.46 

-0.14 

c2dJ033330.4+311051 

3003655 

0.3 

6.6 

0.21 

^•^^-0.26 

62 4+36.di 
0^-‘+_26.81 

2561.06 

339.83 

1 O'*— 
^•^- 0.0 

ii 

- 0.8 

-0.99 

c2dJ034157.4+314837 

3015603 

1.3 

7.69 

0.28 

2.08to 

93.6ll®'^ 

1964.93 

210.13 

^•^- 0.0 

II 

- 0.86 

-1.07 

c2dJ034344.5+314309 

3017083 

2.5 

3.27 

0.28 

3 11+0-3 
^-0.15 

16.84l®®„f 

6740.19 

123.32 

^•^-0.5 

II 

-0.13 

- 0.1 

c2dJ034413.0+320135 

3017363 

4.0 

5.47 

0.23 

3 /|/|+0-3t) 

42.96l?i+i 

17530.01 

245.48 

^•^-0.5 

II 

0.71 

0.61 

c2dJ034441.2+321010 

3007882 

1.2 

4.09 

0.3 

A 0+2-67 
^•^-0.63 

o-43l^:°« 

3965.79 

3561.31 

1 ^,+397.9 
^•■^-0.3 

I 

-0.13 

0.14 


H band turnoff 


c2dJ032519.5-1-303424 

3005058 

5.5 

0.74 

0.27 

1 -T'-j+u.ud 
^•'^-0.06 

40 9 +^^' ^^ 

+^•^-14.26 

311.13 

4.52 

1 A+d.U 

1 -^- 0.0 

II 

0.04 

0.16 

c2dJ033037.0-i-303128 

3000055 

1.7 

8.57 

0.26 

1 4Q+0-^ 
t-+y_0.75 

26 92+00-88 
zo.yz_i3 42 

3108.59 

95.56 

1 1 + 0 -^ 

^ + - 0.1 

III 

-0.76 

-0.92 

c2dJ033044.0-i-303247 

3015197 

1.4 

3.81 

0.24 

A 0+2.34 
^•°-0.47 

. ..+84.64 
°+l^-3.64 

142.18 

145.99 

1 9+5.6 
^+- 0.2 

II 

-0.58 

-0.46 

c2dJ033118.3-1-304940 

3002814 

2.0 

5.21 

0.25 


ji.‘+o_24 58 

202.82 

89.27 

1 c+3.6 
^•^-0.5 

II 

-0.63 

-0.65 

c2dJ033312.8-i-312124 

3011659 

3.8 

8.95 

0.28 

1 .3811^2 

'll '^'1+40.29 

54 

92.66 

507.28 

1 A+0.0 

I 

0.28 

0.03 

c2dJ033341.3-i-311341 

3003019 

2.6 

2.59 

0.23 

l-oo_1.33 

13.841^6,1® 

888.1 

12.61 

1 7+5.1 
^ + -0.7 

II 

-0.36 

-0.29 

c2dJ034109.1-i-314438 

3010092 

3.5 

5.51 

0.26 

2.15t°j« 

49.38lg;^ 

1686.36 

177.4 

1 4+^-^ 
^•+-0.4 

II 

-0.51 

-0.53 

c2dJ034255.9-i-315842 

3019328 

4.1 

5.23 

0.28 

2.29toi^ 

65.84i^°of 

2437.85 

62.43 

c 9 + 0.8 

9 . 9_4 3 

III 

-0.85 

-0.65 

c2dJ034426.7-i-320820 

3011061 

0.6 

0.84 

0.3 

A 21 +014 

n.zi_Q Qg 

l-O^_0.14 

37.06 

17.83 

1 1 + 0.3 
^ + - 0.1 

I 

-0.44 

0.22 

c2dJ034431.1-1-321848 

3005792 

1.2 

3.05 

0.26 

^•^^-0.03 

59 9+38.47 
^^•^-57.78 

3.09 

2.97 

9 7+11.7 
'-2.7 

II 

-0.61 

-0.26 

c2dJ034437.9-i-320804 

3002764 

0.9 

3.92 

0.26 

f\ 0+1.99 
^•^-0.53 

19.6616^1 

93.56 

37.84 

1 9+3.3 
^•■^-0.3 

II 

- 1.12 

-0.95 

c2dJ034516.3-i-320620 

3011667 

2.6 

5.22 

0.23 

1 08 + 1 '^^ 
1.U6_o.28 

71-65ii?+ 

56.7 

49.01 

1 4 + 8-6 
^•+-0.4 

II 

-0.44 

-0.58 

c2dJ034520.5-i-320634 

3011715 

0.4 

4.48 

0.28 

1 co + l-O 

97 09 + 66.1 
Z/. 9 Z 2 i 74 

266.2 

26.14 

1 2 + 2-1 
^+- 0.2 

II 

-0.72 

-0.31 


Ks band turnoff 


c2dJ032741.5-i-302017 

3001576 

0.2 

2.27 

0.24 

0.86+“;^' 

91 99+*^^--*- 
ji.dZ_ 2 i 27 

49.74 

8.74 

1 7+c5.b 

'-0.7 

II 

-0.67 

-0.24 

c2dJ032800.1-1-300847 

3000691 

2.0 

2.53 

0.22 

0-1+0.08 

r ^ + 18.02 

^•^-3.45 

25.08 

7.42 

^•^-0.2 

II 

-0.64 

-0.14 

c2dJ032847.6-i-312406 

3000851 

0.2 

8.46 

0.23 


1 9 qq+63.88 
^^■^^-10.93 

19.85 

1.85 

i-3+o;3 

II 

-1.07 

-0.58 

c2dJ032850.6-i-304245 

3014369 

1.5 

5.59 

0.25 

^+'-0.46 

11.45t®+/ 

36.24 

5.12 

3.4^3+ 

II 

0.04 

0.16 

c2dJ032851.0-1-311818 

3017528 

3.1 

5.9 

0.24 

^•'^-0.96 

45.46t5°;44 

526.49 

66.46 

2+5:o 

II 

-0.65 

-0.34 

c2dJ032851.2-1-311955 

3013451 

1.5 

4.15 

0.23 

1 11+1-06 
^•^^-0.61 

99 1 c+61.02 
jz,.ij_25 55 

169.12 

18.34 


III 

-0.61 

-0.19 

c2dJ032859.6-i-312147 

3015522 

3.6 

2.25 

0.2 

9 9a,+0-82 
^^•^^- 0.88 

13 . 53 I 3 +® 

904.71 

12.41 

^•°-0.6 

II 

-0.77 

-0.55 

c2dJ032903.8-i-311604 

3005745 

2.3 

4.9 

0.25 

0-8+0.65 

^•^'-0.06 

1419.57 

2497.43 

9 7+5.2 
^■'-2.7 

I 

0.99 

1.13 

c2dJ032903.9-i-305630 

3005573 

1.7 

0.65 

0.25 

oi4ig:g^ 

1 19 + 5.31 
^•^^- 1.02 

14.67 

5.3 

19.0_^2.6 

I 

0.01 

0.75 

c2dJ032903.9-i-312149 

3001003 

1.4 

3.29 

0.23 

A QO+l-Ol 

^•^^-0.55 

7 m+69.25 
'•^^-5.74 

206.45 

25.3 

2 

^•^- 1.0 

II 

-0.63 

-0.23 

c2dJ032909.0-i-312624 

3013017 

2.0 

2.11 

0.22 

^+^-0.48 


102.75 

0.72 

9 9 + 5.4 
^•^- 1.2 

II 

0.1 

0.5 

c2dJ032910.8-i-311643 

3005025 

2.7 

1.25 

0.24 

0 3 + 0-21 


11.28 

2.45 

5 2+'^-3 

I 

-0.05 

0.3 

c2dJ032913.1-i-312253 

3012963 

0.2 

2.57 

0.22 

A 27+1-81 

A Q+d.yb 
^•^-0.61 

90.31 

53.94 

^•^- 0.2 

I 

- 0.8 

-0.3 

c2dJ032921.9-i-311536 

3018372 

6.0 

0.77 

0.22 

A 9+0.07 
^•^-0.07 

^•^^-3.01 

34.96 

14.09 

i- 2 ir 2 

II 

- 1.11 

-0.82 

c2dJ032923.2-i-312030 

3009928 

0.3 

0.89 

0.24 

^•^^-0.03 

14 37+27.48 

l‘+.9/_ii 44 

7.74 

2.33 

i+oi 

II 

-0.69 

0.03 

c2dJ032932.6-i-312437 

3003831 

0.8 

1.99 

0.23 

o»-o^_0.53 

13.761^1,5'* 

139.55 

2.06 

l+aa 

II 

-0.91 

-0.38 

c2dJ033001.9-1-303529 

3010577 

4.2 

0.78 

0.24 

0 39+0-13 

0 12+0-31 
'-0.14 

79.42 

1.64 

1 - 8+;8 

II 

0.82 

1.06 

c2dJ033035.9-i-303024 

3009292 

0.4 

3.7 

0.22 

9 'yn+Q.67 
^•^•^-0.44 

39 41+45.77 

1083.2 

125.66 

1 +oi 

II 

-0.94 

-0.87 

c2dJ033038.2-i-303212 

3002146 

1.6 

0.9 

0.24 


0 

ll-4iO_o.i6 

179.08 

170.26 

55.61^^6® 

I 

0.67 

1.76 

c2dJ033052.5-i-305418 

3010241 

0.6 

6.53 

0.23 

-0.46 

91 09 + 6 I .34 

ji.0J_27 74 

44.73 

10.54 
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The YSOs start exhibiting significant (3cr) IR excesses above the photosphere level longwai'd of the turnoff wavebands. 
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0.7 

1.75 
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0.1 

4.07 
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c2dJ034457.9+320402 

3010073 

0.6 

1.86 
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0.04 
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^•^^-0.14 
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Notes: (1): c2d ID; 

(2) : ID of the Best-fit model from Robitaille et al. (2006); 

(3) : Minimum reduced 

(4) : The most probable V-band interstellar extinction; 

(5) : The most probable heliocentric distance; 

( 6 ) : The most probalbe stellar mass and the 95% confidence intei'val; 

(7) : The most probable stellar age and the 95% confidence interval; 

( 8 ) : The best-fit stellar bolometric luminosity; 

(9) : Dust luminosity from integral of the model SED of disk+envelope that best fits the observed SED; 

(10) : The disk inner radius, in units of the dust sublimation radius, and the 95% confidence interval; 

(11) : Evolutionaiy stage; 

(12) : Spectral indices at A > Atumoff for stellar photosphere-included IR SEDs. 

(13) : Spectral indices at A > Atumoff for stellar photosphere-subtracted IR SEDs. 
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